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QQ , 1 Introduction 

OO ' 

' The Newtonian 7i-body problem concerns the motion of n particles with 

, masses G M"*" and positions Qi G M.'^, where i = 1, . . . ,n. The motion is 

' governed by Newton's law of motion 

o^ : 

O ■ qi+li=0 (1) 



> 



where 



^ ■ 7i = ^mkRtkiqi - Qk-) (2) 

cz ' 

and Rik = \\qi — qkW^^. The center of mass qg of the system is defined by 
the relation "^^Ift ~ ^g) = 
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The n particles form a central configuration if there exists a A G R such 
that 7i = X{q, - qa) (Wintner 1941). 

Several aspects of the n-body problem motivate the study of central con- 
figurations. First of all, for suitable initial velocities, central configurations 
may define a homothetic motion (i.e. a motion where only the size of the 
configuration can change) or, in more generality, a homographic motion (a 
motion where the configuration can rotate and change size). In the case of 
a strictly three-dimensional motion homographic motions that are not ho- 
mothetic are forbidden. In the planar n-body problem, instead, there exist 
homographic motions which are relative equilibria (i.e. motions where the 
configuration rotates but does not change size). 

Secondly, for the n-body problem, every motion starting or ending in total 
collision is asymptotic to an homothetic motion and has a central configura- 
tion as normalized limit (Saari 1980). 

Finally, central configurations appear as a key point when studying the 
topological changes of the integral manifolds. The integral manifold I he of 
the n— body problem is the set of points having energy h and angular mo- 
mentum c. Changes in the topology of the integral manifold are caused by 
central configurations. Topological informations on the integral manifold may 
provide qualitative informations about the motion of the bodies. For further 
details and references see Albouy (1993). 

In this paper we will exhibit a new family of spatial central configura- 
tions in the seven body problem. This family is unusual in that it is a simplex 
stacked central configuration, namely a configuration formed by an equilateral 
triangle contained in a tetrahedron (see figure [T]). Examples of trivial stacked 
central configurations (i.e. configurations where a subset of the points is a 
central configuration) are well known. For instance the configurations where 
N particles lie at the vertices of a regular A^— gon are stacked central config- 
urations when TV is not a prime number. Other simple examples are given by 
the "rosette central configurations" , i.e. planar central configurations where 
N particles of mass mi lie at the vertices of a A— gon, N particles of masses 
TO2 fie at the vertices of another 7V-gon rotated an angle of 7r/n from the 
other, and a particle of mass toq lies at the center of the two A— gons (see 
Sckiguchi (2004) and Lei and Santoprete (2006)). 

The first nontrivial example of stacked central configuration, where two 
bodies on a line are contained in an equilateral triangle, was found only 
recently (see Hampton (2005)). The new family of central configurations 
discussed in the present paper seems to be the natural generalization of the 
five body one found in Hampton (2005). 

The results of this paper and those of Hampton (2005) suggest that 
such families of stacked central configurations might generalize to higher 
dimensions, with n + m + 2 bodies in M" arranged as concentric n- and Tri- 
dimensional simplices (for m < n). However, we have been unable to find a 
proof that would easily generalize. 
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2 Central Configurations in M 

We now want to express the equations for central configurations in a way 
that is more convenient for our purposes. From equation ([2]) it is easy to 
obtain 

li - 7i (™i + mj)Rij{qi - qj) + ^ rukiRikiqi - Qk) - Rjkiqj - qk))- (3) 

k^ij 

On the other hand the equation for the central configuration can be expressed 
as 7i — 7j = X{qi — qj). Taking the wedge product of equation ([3]) with the 
vector Qi — Qj we get 

mk{R^k - Rjk)A^jk ^ (4) 

where Aijk = {qi — qj) A (qi — qk)- For a non-collinear configuration, the sys- 
tem of equations (|4]) says that 7^ — jj and Qi — qj are linearly dependent and 
therefore is equivalent to the definition of central configuration. The equa- 
tions ([ll, in the particular case of a planar central configuration, are known 
as the Laura- Andoycr equations and the bivcctor Aijk is simply twice the ori- 
ented area of the triangle {qt, qj,qk)- G. Meyer (Meyer 1933) generalized the 
Laura- Andoyer equations to to higher dimensions (see also Hagihara 1970; 
Albouy 2003). 

Equation ([4]) can be written in a slightly different way if one takes the 
wedge product of equation ([2]), (with the indexes renamed: i = ii, j = 12) 
with ((7,2 — (7^3) A ((7,3 — q,,,) A ... A ((7i^_i — qi^). Thus the equations for the 
central configurations in M"* can be written as 

nik{Rnk - Ri2k)AiA2,....id.k = (5) 

k^ii,...,id 

where = det{qi^ -qi^,... ,qi^_^-qi^,qi^-qk) i.e. it is d\ times the 

signed volume of the d-dimensional simplex formed by the masses ^ , rrii^ , 
■ ■ ■,mi^,mk. 

In particular, when d = 3, the equations above can be written as 

fijh = ^ rnk{Rtk - Rjk)Ajhk = (6) 

where Aijhk is the coefficient of {qi — qj)A{qj — qh) ■ {qn — qk), i-C. it is six times 
the signed volume of the tetrahedron formed by mi,mj,mh-,mk. Moreover, 
since fijh = —fjih, the system of equations ([6]) provides n{n — l){n — 2)/2 
equations. Usually the Aijhk will be replaced by the non-negative volumes 
D'ijhk in order to make the sign of the terms in our equations more apparent. 

Remark 1 If n = c? + 1, and all the mutual distances Vij are equal, that is, 
when all masses form a regular simplex of dimension d, then condition ([4|) is 
verified. Thus the regular simplex of dimension d is a central configuration 
in for any value of the masses. 
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3 Stacked Central Configurations in 

We now investigate configurations with tlic following symmetries and rela- 
tions: 



ri2 


=''13 


= ^23 


= ri4 


= r2i 


= ''34 


ri5 


=?'26 


= ''37 








r45 


=?'46 


= ''47 








r56 


='"67 


= ''57 








ri6 


=''17 


= ''25 


= r27 


= r35 


= ''36 



Figure [T] depicts a typical configuration. These symmetries result in many 
of the Dijki being equal. These equalities will be exploited to simplify a given 
Dijki to Dranop whcrc {mnop} = {ijkV\ and m < n < o < p ; e.g. I?i475 will 
be replaced by 

We want to prove the following 

Theorem 1 For every ris <£ (0, V6/4) there are positive masses mi ~ m2 = 
TO3, TO4 and = ttiq = and distances which satisfy equation ^ with 
1 > ri6 > '"45 > J^is such that the seven points with these masses form a 
spatial central configuration. 

In order to prove the theorem we can start observing that, since the 
configuration is highly symmetric, there are only two degrees of freedom, 
which can be parametrized with ris and rig . Let Sij ~ rfj , then the other 
two distances are determined by 

''56 = - rjr, = si6 - si5 

(that indicates that the polygon of vertices {1, 5, 6, 2} is a isosceles trapezoid) 
and g(si5, si6, S45) := 

-3-|-2si5-3Si5 + 4.Si6 + 4si5Si6-4Sig + 2.S45 + 2si5S45+4si6S45-3.S45 ^ 0. (8) 

The last equation follows from imposing that the 4- volume of the pentachoron 
(4-simplex) of vertices {1,2,3,4,5} is zero. This can be done setting the 
following Cayley-Menger determinant (see Sommcrville 1958) to zero: 
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11111 

1 Si2 Si3 Si4 Si5 
1 S21 S23 S24 S25 
1 S3I S32 S34 S35 
1 S41 S42 S43 S45 
1 S51 S52 S53 S54 



(9) 



where the Sij are computed for the above configuration and V is the 4- volume. 

Recall that, where it is convenient, we replace the Aijhk by the non- 
negative volumes Dtjhk- However in order to do that we need Z\i457 < 0, since 
this oriented volume can have different signs for the configuration defined by 
the equations ([7]). This corresponds to requiring that the triangle 567 is inside 
the convex hull of the regular tetrahedron. 

In the following we will look for the configurations satisfying 1 > rig > 
fAh > fi5- Since the triangle 567 collapses to a point at ris = '\/6/4, we 
restrict our attention to ris =e (0, V6/4). 

Let us condense the above constraints on the configuration into a region: 



Definition 1 ^2 C is the set 

{{ri5,ne,rA5) \ ^15 G (0,^6/4), A 1457 < 0, ff(?'i5,fi6,^45) = 0' 
1 > ri6 > r45 > ris} 

Note that f2 can be written as a graph over its projection in the (ri5, rig) 
plane, and that is usually how we will think of it. In particular, r^^ can 
be expressed explicitly as a function of ris and rig. In a slight abuse of 
notation, we may refer to seven-body configurations 'in S7\ by which we 
mean a configuration whose rij satisfy equation ([7|), equation ([8]), and whose 
{r 15, Tie, r 45) are in f2. 

By assumption we will have ri5 7^ rig, ris 7^ riy. From equation ([5]), with 
z = 1, j = 2, /i = 3, we obtain 

/123 = -m^iRiz - i?25)^1235 - "1-6(^16 " -R26)£'l236 = (10) 

and since i?i5 = i?26, ^16 = -R25 and 151235 = -D1236 we have 

(ms - m6)(i?i5 - R25)Dr235 = (11) 

and we find 7715 = mg. Similarly with i = 2, j = 'i, h ~ I wc find itlq = mj. 
Thus TO5 = mg = 7717. 

Also from equation ([6]), with i = 5, j = 6 and h = 7, wc obtain 

/567 = -mi(i?5i - i?gi)L»i5g7 - 77T.2(i?52 - R(i2)D25(>7 = (12) 

and since D2567 = ^1567 we have 

{mi - TO2)(-R51 - -R6l)-Dl567 = (13) 

and thus mi = 7712. Similarly with i = 6, j = 7, ft- = 5 we find 7712 = nny,. Thus 
77ii = m2 = 7/13. 
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Now let i = 1, j — 4 and h = 7. Then from equation ^ we obtain 

/l47 = m5((i?i5 - i?45)^1457 + (i?16 " i?45)^1467) = (14) 

As long as TO5 ^ the previous equation gives a relatively simple constraint 
on the geometry of our configuration: 

H := (-R15 — i?45)£'l456 — (-R45 ~ ^16)^1467 = 0. (15) 

Thus we are interested in the subset of f2 where H = 0: 

Definition 2 Let Hh E denote the subset of J? for which H = Q. 

Before analyzing H in more detail we prove two technical lemmas. Let 
1/a be the distance between the plane P defined by the points 1, 2 and 3 
and the plane P' defined by the points 5, 6 and 7, and let b = r^Q. With this 
notation we can prove the following 

Lemma 1 Given a configuration in fl we have /147 = 6/241, /175 — ^6/251, 

/645 = 6/461 and /275 6(/251 + /27l)- 



Proof To prove that /147 = —6/241 note that 

/241 = rn^KRie — i?45)(Z\i245 + ^1247) + (^15 — ^45)^1245 

/l47 = — "^5[(^16 — i?45)^1467 + (^15 " -^45)^1457]- 

A simple computation shows that 



(16) 



V36(3 ~V6a + V6ab) 

^1457 = = -OZ\i245 (17) 

_ V36(-6 + V6a6 + 2V6a) _ 

^1467 — T7, — —0(^1245 + ^^1247]- Uoj 

18a 

This concludes the first part of the proof. 
To show that /175 = —6/251 consider 



(19) 



/25I = - "^l(l - -Rl6)^1235 - "14(1 - i?45)^1245 

+ m5(i?i6 - i?56)/\l257 
/l75 = - "li(l - i?i6)Z\l257 + '714(1 - i?45)^1456 

+ "l5(i?i6 — i?56)^1567; 

then it is easy to find the following relationships between volumes 

/il257 - V36/(2a) = -6Z\i235 

\/36(3 - \/6a + \/6a6) .„p.N 

^1456 = 73 = 0ZI1245 y-^^) 

18a 

Ai567 = -V3b^/(2a) = -6Z\i257. 



7 



(21) 



To show that /e45 = 6/45 1 consider 

/461 ="11 [(1 - i?i5)Z\i245 - (1 - i?16)^1247] 

+ TO5(-R45 — -R56)( — ^1456 + ^146?) 
/645 = - {-mi[(l - i?i5)Z\i456 + (1 - i?i6)(/\l456 + ^146?) 

+ mz{RiZ — -^56)^456?} 

then we have the foUowing relationships 

V36(3 -^/Qa + ^/6a6) ^ ^ 

^1456 — 77 — Ozai245 

18a 

x/3b(-3 + ^/6a + 2x/6a5) ^22) 

^1456 + ^1467 = TZ = -0^1247 ^'^'^1 

18a 

_ V3&^(-3 + V6a) _ 

^4567 — — 0(_ZAi45e — Zii467j. 

Da 

FinaUy to prove that /275 = b{f 251 + /271) we consider 

/275 - i?l6)^1257 + (1 - i?l5)-4l257] + "l4(l - i?45)^1467 

+ m5(i?i5 - i?56)Z\i567 

and 

/25I = - mi(l - i?i6)/\i235 - "^4(1 - ^45)^1245 + m^{RlQ ~ R56)^1257 
/271 = - mi(l - i?i5)/il235 - "^4(1 - i?45)^1247 
— m5(_Ri6 + i?l5 — 2_R56)Zli257 

(23) 

Standard computations give the following relationships 

^1257 — — -Ozlli235 

2a 

a/36(6 - 2\/6a - ^/Qab) ^24) 

^1467 = 73 = -0(Zil245 + ^1247) 

loa 

. _ V362 ^ 

^1567 ~ 7 — —0^1257 

2a 

This completes the proof. 

Using the lemma above we have 

Lemma 2 Given a configuration in fi^ there exists positive masses for 
which the configuration is a central configuration. 
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Proof We have seen that, to have /123 = /231 = and /ser — /ers = 0, we 
must choose toi = m2 = TO3 and TO5 = = 7717. Because of the symmetry 
of the configuration and of the above choice of the masses we have that 



/123 = /l24 


= /125 


= /l26 


= /l27 


= /l32 


= /l34 


= /l35 


= /l36 


— /l37 


= /l45 


= /154 


= /23I 


= /234 


= /235 


= /236 


= /237 


= /246 


= /347 


= /264 


= /374 


= /45I 


= /462 


= /473 


= /56I 


= h&2 


= /563 


= /564 (25) 


— h&7 


~ /571 


~ /572 


~ /573 


~ /574 


~ /576 


— /671 


= /672 


= /673 



— /674 — /675 — 

A first inspection of the equations shows that tlie remaining functions that 
must vanish arc /147 and /241, /251 and /175, /271, /46I and /645 and /275- 
However the equations above are not independent for the highly symmetric 
configuration under discussion. Inspecting carefully the remaining equations 
one finds that the remaining functions that must vanish are actually only 
/i47, /251, /271, /461 (see Lemma p. 

Since H = Q, fn-j vanishes. /251, /271 and /461 form a 3 x 3 homogeneous 
linear system of the form Ax = 0, where x = (7711,7714,7715) and in Qu the 
matrix A takes the form 



/ (1 - Rie)Di235 {R45 - l)Dl245 (^^16 - 7?56)I>1257 \ 

(1 -i?15)I>1235 (i?45 - l)i:>1247 (27?56 - i?16 - /?15)/)l257 

V(l -^15)Dl245 +(i?16 - l)I5l247 (J?56 - ^?45) (Dl456 + ^1467) / 

(26) 

For a nonzero mass solution (777,1, 7714, 7775) to exist, the determinant of A 
must be zero. The fact that the determinant is zero can be proved showing 
that the third row is a linear combination of the first two. In order to do 
that, since H = 0, R45 can be written as 



R 



^15^1456 + RieDiiST 



45 



D 



1456 



D 



1467 



so 



(^56 — i?45)(-C'i456 + f 1467) = (^56 " -Rl5)-Dl456 + (^56 " -Rl6)-Dl467- 



Let vi, V2 and V3 be the rows of A. Then a standard computation shows that 
V3 = avi + (3v2 where 



-3 + 2^6 06 + ^6 
9 

3 + -\/6 ab — ^/& a 



■ — —Di24_7 / D1235 
-D1245/ -01235- 



(27) 



This is because 

a{l ~ i?45)^1245 



/3(1 - i?45)i?1247 

-D1247 



-(l-i?45)^1245 



235 



(1 — i?45)£'l247 



245 



(28) 



235 
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Moreover 



since 



(i?56 - Rl5)Di257f3 + (i?56 - Rl(>){Dv257P - I?1257a) 
= (-R56 — ^15)^1456 + (-R56 — Rl6)Dl467 



(29) 



x/36(3 - V6a + x/6a&) _ 

^1257P — 73 — ^1456 (^UJ 

-L'l257P - ^'12570: = j^^^ = -D1467 (31) 

Finally the fact that the masses may be chosen all positive is a consequence 
of the sign pattern of the coefficient matrix. For the configurations in fi^ 
the sign pattern is 







(32) 



where * denotes a sign that is unimportant for the discussion. The first 
and last row imply that all the components of a nonzero nuUvector have the 
same sign. 

To complete the proof of the theorem we only have to show that for every 
ri5 e (0, •\/6/4) we can choose rie and r45 so that H = and (ri5, rig, r45) G 
f2. Note that in f2, H is negative when D1456 = (i.e. Z3i457 ~ 0, since 
-01457 = -Di456) and when R45 = i?i5. H is positive when R45 = Riq. For a 
given value of < ris < r^g one can vary rig from the curve R45 ~ Riq to 
the curve -D1457 — (Figure 2). Thus _ff = is satisfied for some value of 
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rig. For rl^ < ri5 < •\/6/4 (i.e. when the corresponding point on the curve 
-R4.5 — R15 is in Q) one can vary rig from the curve R45 ~ Rie to the curve 
-R45 — i?i5 = 0. Hence H = for some value of rig. 

Therefore there is at least one central configuration for every ?'i5 € 
(0,^6/4). 

We can show a little more: 

Theorem 2 There is only one such configuration for each ri5 G (0, \/6/4). 

Proof To prove that there is only one such configuration for each ris it is 
enough to show that H' < 0, where the prime denotes differentiation with 
respect to rig and with ris held fixed. From equation ([5]), using implicit 
differentiation, we have 



r,,{l + rl,+2rl,~3r%) 



(33) 



which is negative for the configurations in i7. 

Recall that H = {R15 — i?45)£'i456 — {R45 — Riq)Di4qt. Now consider the 
derivative 

H' = 3rj5*r45Di456 + (-R15 - i?45)£''i4.56 + i^f^b'^ib ^ 3rj"g*)Di467 
— (i?45 — -Rl6)-D'i467 

==9n^ ^^?456 + (^15 - i?45)(i??456)') ^^^^ 

^^^1456 

1 -'„-4, 



The point of the above rearrangements is to get all of the volume expres- 
sions into polynomial form. The only parts of this expression whose signs are 
not clear on the region of interest are the volume derivatives (^^1467)' and 
(Di45g)'. Unfortunately, the latter does indeed become negative within Q, 
and even on the set = 0. 

We will first show that (i^ugr)' is positive on the set Hh- Let us write 

/r,2 V _ 4ri6(r^g -r^5)ti 
1^1467 J - 1 , „2 ,02 _ 0^2 ' 

where ti is a polynomial in ris, rig, and r^^. The sign is clearly determined by 
the sign of ti. If we compute the remainder t2 of ti divided by the polynomial 
in equation ([8|), with respect to the lexicographic monomial ordering rig >- 
fib >~ fib, we find that 

h = {l- r?5 + ds) [2(r?6 - rl,){l - rj, + rj,) + UDl^,] , 

where D145 is the area of the triangle formed by the points qi, (74, and (75. 
This form of the expression is clearly positive. 

To rule out a zero of H' on the set H = we began by computing 
two Grobner bases (using Magma) of the following systems. Both were con- 
strained to include the set H — 0, hy using the following two polynomials: 
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H multiplied by {R15 — i?45)£'i456 + {R45 — ^16)^14677 with denominators 
cleared (denoted below as 61), and the condition ([8]) converted to the 
variables (denoted as 62). Both of these polynomials also possess a number 
of components irrelevant to our analysis, but this is unavoidable. To exclude 
components of the solutions of these systems with any of the = 0, the 
condition ri5ri6r45W — 1 = was also added (denoted by 63). These three 
polynomials will be referred to as the base system: 



/i, — ^6 I 3 3 10 I 8 8 6 8 I 2 6 8 a 5 3 8 6 2 8 

Ol — '45' 15 + ^'16' 45' 15 + '16' 15 ^ '16' 15 + ^'16' 45' 15 ~ ^'16' 45' 15 ~ ' 16M5' 15 

_ o^6 3 7 _ 10 6 _ 10 6 I ^8 6 , r, 2 8 6 , r, 8 6 , r, 3 7 6 
^'^16M5' 15 ' 16' 15 M5' 15 + '16' 15 + ^'16' 45' 15 + ^'45' 15 + ^'16' 45' 15 

— 4rfgr|5r45 + '2.r\^T\^T\^ + r\^r\r^r\r^ + r-fgr^grfg + 2r\^r\^r\^ + 2r\^r\r^r\r^ 
+ 2rigr45r^5 + r^grfgrfg — 2r^grj5ri^5 + 2r\(^r\r^r\v^ + 2r^gr45r45 — 2r^^r%^r\r^ 

I 9^.8 T-3 T-3 Ot-S 7.3 ^6 8 2 8 6 2 6 6 2 , 6 10 

"r^'l6'45'l5 ^'16'45'15 'l6'45'l5 'l6'45'l5 ' 16 ' 45' 15 + ' 16' 45 



„8 8 6 8 I 10 6 8 6 , 6 6 

16' 45 ~ ' 16'45 + ' 16'45 ^ '16' 45 + ' 16'45' 



h = n^rieri^w - 1 



45 



(35) 



The first system also included the sign-determining factor of (ir^^r'^c,Dl^;^g+ 
(i?i5 — -R45i(£'i456)'j ^ polynomial pi of total degree 11 in the variables ris, 
rig, and r45E|: 



pi = 6rll - 24r2gr?5 - rfgr^g + SGrf^,, + Sr^^r^g - brl.rj.rl - 16rl,rl 
+ ''45''i5 ~ SOrfgr^g — 3r45rj'5 + ISr^gr^g — ler^grlgr^g — 2r|5ri5 

- ISr^grfs + 22r4gr25rf5 + 22rl,rlrl + r^^rf^ + Grf^ - Srlrf, 

+ Sr^grlgrfg — 2rl^ri^ + 12rfQrl^ + 6r45ri5 — ISrf^rl^ — ISr^grlgrfg 

- 8r|5r?5 + 18r4gr?5 + 26r4gr45r?5 + Sr^rl, - Gr^gr^g - 22r%rl^rl^ 
+ "^rterlrh " Sr^grlgr^g - 2rlrl, + 9rlrl - 2rl,rlrl + 2rlrl, 

- Arf^rl^rf^ — 4rigr|5ri5 — rfgr^g — 6rll + Or^grfg + Srfg — Sr^grjg 

- Srjg + Ar%r% - Artevl^ + Sr^gr^g + 2r% 



Since H' has the form —(a + bpi) with a, b positive functions, if pi is positive 
on S7h then H' < on J?//. After eliminating w, a lexicographic term order 
with r45 ;^ rig ;^ ri5 was used to find a polynomial qiijir-,) of degree 404. 



^ To speed the computation, the actual initial system used in the Groebner basis 
computation included another polynomial in the ideal generated by pi and 62, but 
this is theoretically redundant. 
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The first and last few terms of qi are shown below: 

gi = r^04 _ (804495232776247942366788446548795369/ 
6172585442506994000774739535727460)r^g^ 

+ (19481733600623784042235202979633272487804463/ 
2759886403053727157626401541214461915200)rfg° + . . . 

- (139/99902129761265947794751686948283396618716394291200)ri5+ 
(1/66601419840843965196501124632188931079144262860800) 

(37) 

There are four zeros of qi in the interval (0, ^); a rigorous demonstration 
of this is possible using interval arithmetic. These zeros are approximately 
at n = .5104 ...,r2 = .5384 . . ., = -5774 . . ., and r4 = .5856 . . .. 

For the second system, we first eliminated I?i467 from the expression of 
H' by using the condition that H = rearranged as I?i467 = -Duse^^E^^- 
After factoring out the positive quantity I?i456 we can write 

H' = [6r45V45^L56 + (^15 " ^45)(^?456)' + 

-n.15 — U45 

Now we clear denominators from the following partial sum: 

(i?15 - i?45)(i^?456)' + (^ f'~f' ir^5r'45 ^ -lt)^?467 
iti5 — ^45 

and add the numerator p2 to our base system. If this polynomial does not 
vanish on Hh then neither docs H' . The full polynomial p2 is: 



- '2^r%rl,rll + 60rJ«rl| - SrHrH + 6Qr%rll - ArlrlrH - Ar^rll 
+ Wrf.rlrll + 2Ar%rll + 8rl,rl,rll + 8rl,rl,rll + I6rt,r%rll 

+ 56rl,rl,rll - ir%rlrll + 24r%rl,rll + GOrf.rl.rll + 36r%rl,rll 

■JU'^15'^16 + ^^45^^16 •J^^is'^ie ^^15'^45'^16 J-°U5M5M6 ^Ur^gfig 

4_ c„4 8 10 _ A 2 8 10 , o 8 10 , a 5 7 10 _ o 6 10 , o„6 6 10 
+ ^' 15' 45' 16 ^'15' 45' 16 + 'JM5' 16 + ^' 15M5^16 O' 15' 16 + •J' 15' 45' 16 

-'^'J'^15' 45' 16 + "'15' 45' 16 O' 15' 45' 16 'J^' 15' 45' 16 + "-''15' 45' 16 

- 25rltrlrl', + 8r%rlrl', + 9rlrl,rl° + 6r%rlr% - 42r%r%r% 

+ Gr%rl,r% - 60rllrl,r% - 60rl,rl,r% - 2Ar%rl,r% + 6rllr% ~ Grltrl, 
+ 6rllr% + 9rl,rlM, + 8rllr% + 12rl,rllr% + 6rllr% - 8rt,rllr% 
+ 2ri5r45rig — 8r]gr^g — Wrl^r^^rf^ — IGrf^rlr^rl^ + Grl^rf^ + Trf^rl^rl^ 
~ 2''i5''45^i6 ~ ''i5''45^i6 + 2?'45?'i6 + IS^'is^'Is^'ie — 4rf5rJ5r^g + IQrfr^rl^rlQ 
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^' 15M5'16 ^' 15'45' 16 ^'15'45' 16 "'^ISMS' 16 + 15M5'^16 

14ri5r4gr^g — Trf^rlr^rfg + 2rlr,rlr^rfQ — ISrigrfjr^g — 12ri5rJ5r[g 
6r?5rl5r[g + 24r}°r|5r[g - ^tX^^^^^ - (Sr^r^rl^ + ^^r^^r^rl^ 

6''i5'"45fi6 + 9^i5^45''i6 + l^rfgrfgr^g — 2\r\r^r%^r\Q + Vbr^T\^T\ 



15' 45' 16 "T 15' 45' 16 "r ^^' 15' 45' 16 ^^' 15' 45' 16 "r ^^' 15' 45' 16 

+ 12r^5rj5r^g + ISr^grJgr^g — Or^gr^gr^g — 3rjgr|grfg + Sr^grlgT^g 

- Sr^grlgr^g - ^rWrl^^ - (ST^^r\^r\^ - ^r^^rl^^^ - 6r?gr|gr^g (38) 

+ 2^T%r'gr\^ + 42r8grl0r4g + erJOr^gr^g - Grfgrfgrfg - Qrfgrlirfg 

+ Srl^rigrjg + 21r8gr9gr3g - grj^rjgrfg - ^r^^rl^r\^ - ISr^grlgrJg 

+ 3''i57'45'"i6 + Srjgrlgrjg + 3rJ5r|5r4g + irX^/^r^rX^ — ^r\yj^r\Q 

_01„8 12 2 I 01 „6 12 2 _oi„10 10 2 _ fi„8 10 2 „ 1 c;„6 10 2 
^-'-' 15'45' 16 + ^-'^' 15'45' 16 ^^'l5'45'l6 "'15'45'l6 ^3'l5'45'l6 

+ Srlfrt^r^g + Xhr^^rl^^^ + Zrl^^^^^ + Sr^gr^gr^g + gr^^rl^ 

10^12 01^8 ^12 , 0^12^10 , R^IO^IO , irJ ^10 Q^14 S 



'5'^15'^45 '^'^^15^45 + ■^^15'^45 + "'^15^45 + -Lf^l5'^45 •J^15'^45 

Vl2„8 _ o 10 8 _ o 8 8 
■J'lSMS "J'15'45 •J'15'45 



Again we computed a Grobner basis, in the same manner as for the first 
system, and obtained a polynomial (72 (''15) of degree 466 which is in the ideal 
generated by the polynomials of the second system. The first and last few 
terms of q2 are: 



52 = ? if - (397970325260732999941673276305737323331716006483/ 
4063279255141998390265984630470563259415746780)r^f 

+ (717654067077578250690275945137718952732059392991627077/ 
585112212740447768198301786787761109355867536320000)r^f + . . . 

+ (769475777092319973480299/2698353810706234317867355469145285098 
555299416887370190209448673280000)rl52 

- (2505469532410439724481/8994512702354114392891184897150950328517 
66472295790063403149557760000) 

(39) 

Again it is possible to rigorously compute the existence of exactly three 
roots in the interval (0, -^). These are located at approximately r^, — .5004 . . ., 
re = .5027.. ., and ry = .5252.... 

If we consider the arrangement of the roots of the polynomials qi and 
q2, wc see that we need only verify that on the set J7_y, pi is positive for 
ri5 e (0, ri), p2 is positive for ri5 S (rr, ^), and that either pi or p2 is 
positive for ris e (rijTy). To do this, we can pick three rational values for 
ri5 each of the three intervals just described and compute the Grobner basis 
of the base system specialized to those values. We chose to use the values ^, 
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Fig. 3 Stacked central configuration family 



1115 ~ 1 




(a) (b) 
Fig. 4 (a) Masses of the configurations, (b) Close up of masses near ms = 0. 



ii, and |. Using interval arithmetic, it is possible to verify the appropriate 
signs of pi and p2 ■ 



4 Numerical visualizations 

Since it is often helpful to have an accurate sketch of phenomena, we repro- 
duce here an image of the family of central configurations described above 
(Figure[3|), and the associated masses (Figured]). The masses were normalized 
so that mi + + ms = 1 and then projected into that plane. 
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